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Dynamical screening effects in correlated materials: plasmon satellites and spectral weight 
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Dynamical screening of the Coulomb interactions in correlated electron systems results in a low-energy ef- 
fective problem with a dynamical Hubbard interaction U{uj). We propose a Green's function ansatz for the 
Anderson impurity problem with retarded interactions, in which the Green's function factorizes into a contribu- 
tion stemming from an effective static-U problem and a bosonic high-energy part introducing collective plasmon 
excitations. Our approach relies on the scale separation of the low-energy properties, related to the instantaneous 
static U, from the intermediate to high energy features originating from the retarded part of the interaction. We 
argue that for correlated materials where retarded interactions arise from downfolding higher-energy degrees of 
freedom, the characteristic frequencies are typically in the antiadiabatic regime. In this case, accurate approx- 
imations to the bosonic factor are relatively easy to construct, with the most simple being the boson factor of 
the dynamical atomic limit problem. We benchmark the quality of our method against numerically exact con- 
tinuous time quantum Monte Carlo results for the Anderson-Holstein model both, at half- and quarter-filling. 
Furthermore we study the Mott transition within the Hubbard-Holstein model within extended dynamical mean 
field theory. Finally, we apply our technique to a realistic three-band Hamiltonian for SrVOa. We show that 
our approach reproduces both, the effective mass renormalization and the position of the lower Hubbard band 
by means of a dynamically screened U, previously determined ab-initio within the constrained random phase 
approximation. Our approach could also be used within schemes beyond dynamical mean field theory, opening 
a quite general way of describing satellites and plasmon excitations in correlated materials. 

PACS numbers: 71.27.+a, 71.30.+h, 71.10.Fd 



I. INTRODUCTION 

Over the last years, significant progress has been made in 
the modelization of strongly correlated materials. Such sys- 
tems typically contain partially filled rf or / orbitals'-^, which 
lie relatively close to the nuclei. Electronic Coulomb interac- 
tions can then induce substantial corrections to a Bloch one- 
particle picture, ranging from renormalizations of effective 
parameters in the sense of Landau to full localization of the 
d- or /-degrees of freedom in the Mott insulator. 

One of the difficulties in describing such correlation effects 
is to separate the usually rather small energy range of the cor- 
related d and / orbitals from the larger energy scale of the itin- 
erant degrees of freedom, e.g. the p-orbitals of ligand atoms, 
but also of higher or lower lying states of the transition metal, 
rare earth or actinide atoms themselves. The former orbitals 
are mainly responsible for the low-energy physical properties 
of the compounds, while the latter act as a screening medium, 
setting in particular the actual value of the Coulomb repulsion 
of the correlated degrees of freedom. 

As discussed in Ref. 3, screening is a dynamic process 
which leads in general to a frequency dependent U — U{lo). 
In a realistic approach to strongly correlated materials, U 
can be determined at the random phase approximation (RPA) 
leveli, once the bands and their eigenstates are computed by 
an ab-initio calculation. The static value Uq — U{{)), evalu- 
ated in this way^'^, has been recently used in the dynamical 
mean field theory (DMFT) calculations of materials^, and the 
effect of the frequency dependent screening has been either 



neglected or empirically taken into account by adjusting the 
effective static Uq. The DMFT approach, combined with the 
density functional theory (DFT), is an extremely powerful tool 
to treat ab-initio strongly correlated systems, once the low- 
energy model is determined^. However, very little is known 
on the impact of the frequency dependence of the interaction 
in the low-energy part of the spectrum. The hardest obstacle 
in order to include the dynamic U into the DMFT framework 
has been the lack of a reliable solver for the quantum impurity 
problem with a frequency dependent Hubbard interaction. 

The [/(cj) computed so far for various materials^^ shows 
some common features. The unscreened U (Uoo) is up to an 
order of magnitude larger than the screened one (Uq). The 
frequency dependence, although complicated, can be repre- 
sented roughly by a single plasmon frequency wq, whose 
value is usually much larger than the bandwidth (antiadi- 
abatic plasmon). If one wants to deal with its precise 
form, U{uj) can be resolved in many plasmon contributions, 
which characterize the screening process. These features 
make the problem difficult: The large Uoo rules out the ap- 
plication of traditional weak coupling expansion methods, 
while the presence of many plasmons prevents the direct 
use of methods developed in the context of the Hubbard- 
Holstein Hamiltonian^"'^, as many bosonic baths coupled to 
the fermion degrees of freedom will be necessary to fully re- 
solve U{uj). Recently, this problem has been overcome by 
a continuous time quantum Monte Carlo (CTQMC) solver 
proposed by Werner and Millis'^'^, where a multi-plasmon 
Lang-Firsov transformation''^ is treated exactly in the con- 



text of a hybridization expansion algorithm for the DMFT 
impurity Hamihonian^"'^'. Also the weak coupling CTQMC 
algorithm^^-^ by Rubtsov can treat generic retarded interac- 
tions, but it is limited to a not-so-large dynamic U and not- 
so-large screening frequencies, and therefore it becomes pro- 
hibitively costly for realistic applications of dynamic screen- 
ing interactions in a multi-orbital context. 

Another major problem still left (even if realiable Monte 
Carlo data are available) is the possibility of computing spec- 
tral properties, such as high-energy plasmon satellites. In- 
deed, the presence of screening modes leads usually to a quite 
complicated spectrum with a series of peaks located at high 
energies (at multiples of the plasma frequencies). Those fea- 
tures are hard to get by the usual Maximum Entropy (ME) 
methods^"*-^^, used to invert the noisy QMC data in the imagi- 
nary time domain into the spectral properties at real frequen- 
cies. The ME methods are quite reliable at low energy, but 
usually are not capable to deal with high frequency features. 

In this paper, we present a DMFT approach based on a Bose 
factor ansatz (BFA) for the Green's function which is able to 
handle a generic U{uj) interaction in a strong coupling antia- 
diabatic regime, a typical situation in strongly correlated ma- 
terials, and provides a robust and general way to compute the 
full spectrum of the frequency dependent (retarded) U, with 
an accuracy capable to resolve the high energy satellites. Our 
method is based on the separation between the energy scales 
set by the screened value Uq and treated using well established 
solversii22i2^, and the dynamic part treated with various levels 
of approximation, the simplest and most insightful one taken 
from the dynamic atomic limit. 



The paper is organized as follows. In Sec. |ll] we spec- 
ify the Anderson impurity model we would like to solve and 
the Green's function ansatz used in our method, in Sec. |lll| 
we present the dynamic atomic limit approximation (DALA) 
to our approach, in Sec. |IV] we show various ways to im- 
prove upon the DALA, and in Sec.[V]we describe their perfor- 
mances. In Secs. lVII and lVIII we report our results for a single- 
band lattice model at half-filling and a three-band model with 
the DFT density of state (DOS) of SrVOa, respectively. Fi- 
nally, Sec. I Villi summarizes our findings. 



II. THE GREEN'S FUNCTION BOSE FACTOR ANSATZ 
(BFA) 

A. General model 

We discuss here the case of a multi-orbital Hubbard- 
Holstein model. In the context of the DMFT approach, one 
maps the full lattice problem into a single-site Anderson im- 
purity problem coupled to an effective bath. The bath is deter- 
mined self-consistently by requiring that the impurity Green's 
function equals the on-site Green's function on the lattice^i 
Therefore, computing in the most effective way the Green's 
function of the Anderson model (AM) is of key importance 
to have a feasible DMFT scheme. Here, we have the addi- 
tional complication that the effect of screening makes the on- 
site Hubbard interaction retarded. In the Matsubara imaginary 
time action formalism the dynamic Anderson model reads: 



J 



S 



jy jyY^clir)g,l{r-r')cUr') + l J^^dr jy ^ (^n,.(r) - 1) C/,.,,..(r-r') (^n,,,{r') - 1 



(1) 



where /3 is the inverse temperature, Uia-.ja' (t) is the screened 
interaction, {i,a) is the set of orbital and spin indexes, and 
(c|^, Cia) are the creation and annihilation operators satisfy- 
ing the antisymmetric commutation relations. QQl^{iu}) — 
iid + II ~ Aia-{iuj) is the effective hopping term coupled with 
the bath via the hybridization A,;(j(Jw), and /i is the chemical 
potential. Let us assume that the dynamic part of U is orbital- 
independent and couples only to the total charge. 



C/,.,,<x'(r)=f/P ,<,,<S(r) + C7(r), 



(2) 



with U{0) = 0. Therefore, we assume that all details of the 
interaction due to effective spin and Hund's couplings are em- 
bedded in the instantaneous part of the full U. Here and there- 
after we are going to take the convention that the instanta- 
neous t/'^ is the static (fully screened) limit of U with U being 
a repulsive contribution. 



B. The Bose factor ansatz 

We are mainly interested in evaluating the Green's function 

Gij{T) — (Tcj(r)c|^(0)) and its spectral properties for the 
model in Eqs.[T]and|2l We are going to rewrite it in the form: 



G,,(t)=F,,{t)Go^,,(t), 



(3) 



where Go,ij (r) is the Green's function for the model in Eq.[T] 
but with a static on-site repulsion, namely Uia-,ja'{T) — 
^ia ja-'^ {''')■ ^^ highlight that the above factorization is de- 
fined in the time domain, a feature which is borrowed from the 
dynamic atomic solution of the problem whose form is known 
analytically, as explained in Sec. |III] In that limit the Green's 
function assumes exactly the form in Eq. [3] with Go the in- 
stantaneous Uq atomic Green's function. The static model is 
much easier than the dynamic one, since it contains only the 
energy scales set by the screened U{-^ Uoo) and the Kondo 
resonance with the bath, and it can be solved by means of 



various techniques 



17.22.26 



which are usually very robust and 



efficient in this case. On the other hand, Fij (r) is a Bose 
factor, which is a functional of U{t), and it is not known 
a priori. However, we will present various approximations 
where the function Fij (r) is derived. It contains the informa- 
tion of the plasmon (or phonons) excitations, and the plasmon 
(or phonons) satellites. 



C. BFA spectral properties 

A great advantage of deaUng with the Green's function 
ansatz in Eq. [3] is the possibility to compute very accurate 
spectral functions over the whole energy range, including the 
intermediate-high energy plasmon satellites. Indeed, since the 
Bose factor Fij (t) can be estimated analytically by means of 
some approximation, its numerical value is known at machine 
precision, and its spectral function B{uj) can be obtained via 
a Pade approximant,^^ in an accurate and robust way. On the 
other hand, a ME approach has to be used to find the spectral 
function Ao{uj) of the static Green's function Go,ij{^)- How- 
ever, this does not pose any particular problem, since there 
are no high-energy features in Ao{lu), and its energy range 
is set by Ua, where the ME is reliable in presence of data 
with good statistics. ^^ The spectral function A{ll!) of the full 
Green's function G^ (r) expressed as a functional of B and 
^0 reads: 



A{io) = / de B{e) 



1 + e-^" 



{1 + e-f>(^~'^^){l - e~P^) 



Ao{uj~e). 



(4) 

The spectral functions obtained in this way are reported for 
instance in Fig. |2] which corresponds to the Green's func- 
tions plotted in Fig.[T] The quality of the satelhte resolution is 
striking, much higher than the one usually obtained with ME 
methods^'', particularly at energies far away from the Fermi 
level. 

We would like to stress that the spectral convolution in Eq.|4] 
is general, and can be used not only for the approximated 
Green's functions we are going to derive in Sees. Hill and IIVI 
For instance, our approach to compute the spectral properties 
can be applied to the Green's function obtained by means of 
the algorithm in Ref. 2 1 . Given the full Green's function G{t) 
of the dynamic impurity problem of Eq. [T] one defines an 
auxiliary Green's function GauxlT") as G{t)/Fuala{t), with 
^DALA taken from the atomic limit as described in Sec. [Ill] 
This is an effective way to exploit the separation of the low 
energy properties, kept in Gaux, from the high frequency fea- 
tures correctly reproduced by the DALA. At this point, one 
computesylo, the spectral representation of Gaux, by using 
ME, and B, the spectrum of i^DALA, by means of the Pade 
approximant, and evaluates the full spectral function in Eq. |4] 



in. THE DYNAMIC ATOMIC LIMIT APPROXIMATION 

In order to find a way to determine Fij (r) in Eq. [3j for the 
moment let us take into account the single-orbital symmetric 



case. Thus, we can drop all the orbital and spin indexes, and 
simplify considerably the notation. By inverting Eq. [3] one 

gets: 



F(r) 



Gir) 



Go{t) \Go{r) 



G{r) 



(5) 



A=0 



where the rightmost-hand side of the above Equation is the 
approximation for F{t) taken in the dynamic atomic limit 
(DAL), when the hybridization A is zero. It turns out that 
the Green's function in the DAL is analytically solvable by 
means of a Hubbard-Stratonovich transformation^**, and there- 
fore G/Go is exactly known in a close analytic form, such 
that: 



-FbALA(r) = exp 



E 



U{iVn) - Uq 



1) , (6) 



where Vn ~ 2n7r//3 are bosonic Matsubara frequencies, with 
n relative integer 

Beside the atomic limit, this approximation is exact in 
the static and the non-interacting limits (in both cases 
Fbala{t) — 1). Notice that it retains all the non-perturbative 
character of Go(t). To have a better idea on the quality 
of this approximation, we are going to test it for the dy- 
namic U with a single plasmon mode wq, which is equiva- 
lent to the Anderson-Holstein model with 'U{iv) — Uoo — 
2X^100/ {ly'^ + ojq) and the electron-"phonon" coupling given 
by A = y^{Uoo~ Uo)loo/2. We use the CTQMC algorithm 
by Rubtsovr^i^^ which can handle retarded interactions and 
yield the exact Green's function in a weak coupling regime, to 
benchmark our approximation for the particle-hole symmetric 
system with Uq = 2, /3 ~ 10, and few values of wq and Uoo- 
The energy units are expressed in terms of the half bandwidth 
(D/2 = 1) of the semicircular DOS. 

As one can see from Fig. [T] the DALA works very well 
for the cases analyzed, since it gives a Gdala(t) which al- 
most coincides with the numerically exact G{t) given by the 
CTQMC algorithm. The accuracy is particularly impressive 
in the case f3 = 10, Uq — 2, Uoo = 6, and cjq — 5, where the 
impact of the dynamic part is reduced by the larger wq ( 3> Uq), 
and a smaller Uoo, namely when the the energy scales of the 
static part set by Uq are well separated form the dynamic con- 
tributions in U{iv). Moreover, from Fig.[T]it is apparent that 
the low energy properties of the system are strongly renormal- 
ized by the effect of the high energy components of U . This 
is quite insightful on the importance of the dynamic screen- 
ing effects in the treatment of more realistic models that we 
are going to tackle in Sec. |Vll This can be noted also from 
the spectral properties reported in Fig.|2] which correspond to 
the Green's functions Gdala(''') plotted in Fig.[T] A spectral 
weight transfer from the low frequency spectrum to the high 
energy satellites is clearly visible in the Figure. 

Since the DALA is obtained in the A = limit, it works 
well in the intermediate-strong coupling regime, with C/o and 
the dynamic part large, as we have seen in the cases ana- 
lyzed in Fig. [T] where Uq was quite close to the critical Uc2 
(sa 2.6, see Ref.ll^) for the Mott transition of the static Hub- 
bard model. However, it deteriorates as Uq is getting smaller 
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FIG. 1: (color online) Green's functions for a half-filled Anderson 
model with dynamically screened U and a semicircular density of 
states at /3 = 10. The Gdala obtained by the method proposed in 
Sec.|lll]is plotted (green long-dashed line) and compared to the exact 
Rubtsov's CTQMC numerical result (blu dotted line). Also the dy- 
namic atomic limit (red solid line) and the static one (pink dot-dashed 
line) are reported. The system is computed at a quite large tempera- 
ture, such that an accurate benchmark against the numerically exact 
CTQMC is still possible, even for quite large values of U. 
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FIG. 2: (color online) Spectral function of the Gdala Green's func- 
tions reported in Fig.[T] 



and the hybridization A becomes important to set the low- 
energy properties of the system. To show this, let us take into 
account the Anderson-Holstein model with Uq = 2, loq — 5, 
and Uoo — 6 at /3 — 10, for which the DALA gives a re- 
suh very close to the exact one. Now, let us keep loq and 
Uoo — Uq fixed, such that the DALA Bose factor (Eq. |6]l is 
unchanged, while we vary Uq from strongly to more weakly 
correlated values. The resulting Green's functions are plotted 
in Fig. [3] It is apparent that the DALA accuracy reduces as 
C/o decreases, and the resulting bias is more pronounced in 
the low-energy part of the Green's function. Indeed, the dis- 




FIG. 3: (color online) Half-filled Anderson impurity model with re- 
tarded screened interaction: uo — 5, Uoo — C/o = 5, /3 = 10, for 
various values of Uo (a) numerically exact Rubtsov's CTQMC and 
DALA Gir); (b) CTQMC and DALA imaginary part of G(itJ„). 



crepancy is larger around [3/2 in G{t) (Fig.[3t), which corre- 
sponds to a larger difference at small Matsubara frequencies 
in G{iujn) (Fig.IlJ)). On the other hand, the high-energy tails 
of G{t) are very well reproduced by the DALA, as it is con- 
firmed also by the inspection of G{iuj) at large a;„(> 4). 

The correct high-energy asymptotics of the DALA is a non- 
trivial property of this approximation, which is borrowed from 
the atomic limit exactly built-in. In order to further analyze 
this important feature, we take into account the temperature 
dependence of the DALA in the symmetric Anderson model 
with a dynamic interaction given by Uq = 1.25, ujq = 2, and 
Uoo — 2.5. In this not-so-correlated case, the model can be 
solved exactly down to low temperatures (/3 = 160) even by 
the CTQMC algorithm, to benchmark the temperature depen- 
dence of our approximation. Results are plotted in Fig. 2] As 
we already found in the previous analysis, at large Matsub- 
ara frequencies the dependence of the Green's function is cor- 
rectly given by the DALA, which in this case becomes almost 
indistinguishable from the exact CTQMC result for a;„ > 7. 
We note that the DALA is capable to reproduce the decay of 
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FIG. 4: (color online) Symmetric Anderson model with retarded 
screened interaction: Uo = 1.25, a;o — 2, f/oo ~ 2.5, at different 
/3. (a) Imaginary part of the numerically exact Rubtsov's CTQMC 
and DALA Green's functions for small Matsubara frequencies uun', 
(b) —ijJnlm[G{iujn)] in the intermediate frequency range. The arrow 
indicates the frequency when the DALA and the CTQMC Green's 
functions become practically indistinguishable. The DALA repro- 
duces correctly the intermediate-high energy behavior of the Green's 
function well beyond the "trivial" l/iuj term; (c) G{t); (d) The 
CTQMC and DALA Bose factor F{t). The thicker lines are for the 
DALA. The discrepancy between the DALA and the exact factors is 
getting smaller as the temperature increases. The wiggles are due to 
the stochastic noise of the data. 



the imaginary part of G(ici;) well beyond the l/ioiterm, as it is 
apparent from Fig.|4fb). It is also worth noting that the relative 
accuracy of the approximation increases with the temperature, 
as it is shown in Fig.Ud). Indeed, the exact Bose factor F{t) 
is getting closer to Fdala(t) as the temperature increases. At 
low temperatures, it is the F(/3/2) value which is poorly re- 
produced by the DALA. Again, this is related to the rough- 
ness of the approximation at low-frequency, which does not 
describe accurately the low-energy excitations around and be- 
low the coherent temperature. Indeed, the Friedel sum rule is 
clearly violated, as one can see in Fig.Ua), where the condi- 
tion G (iO+) = ~4/D valid at half filling is not met by the 
DALA. 

Therefore, going beyond the DALA is needed to capture the 
low-energy low-temperature features of the spectral function, 
while its high-energy properties, as the plasmon satellites, can 
be successfully taken into account at this level of approxima- 
tion. 



IV. BEYOND THE DYNAMIC ATOMIC LIMIT 

A. The DALA Lang-Flrsov approximation 

A way to improve the dynamic atomic limit in both the 
intermediate-low energy correlations and low temperature 
regimes is provided by the Lang-Firsov approach. This ap- 
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FIG. 5: (color online) Symmetric Anderson model with retarded 
screened interaction as in Fig. |4l Uo — 1.25, ojo = 2, f/oo = 2.5. 
(a) Imaginary part of the numerically exact Rubtsov's CTQMC and 
DALA-l-LF Green's functions for small Matsubara frequencies (jj„ 
at different 13; (b) The CTQMC, DALA and DALA+LF Bose fac- 
tor F{t) at /3 = 160. The improvement of the F{t) provided by 
the Lang Firsov correction is apparent. The wiggles are due to the 
interpolation of noisy QMC data; (c) Spectral representation of the 
DALA-l-LF Green's function in (a). 



proximation has been widely used in the literature to tackle 
electron-phonon models in the antiadiabatic limit, when the 
electron-phonon coupling is A ^ uq, with ljq the phonon fre- 
quency. The same applies to models where the interaction is 
retarded by the charge screening plasmons. In the latter case, 
the antiadiabatic regime is met more often, since the plasmon 
frequences are larger than the phonon ones, usually by an or- 
der of magnitude. However, it should be noted that in the real- 
istic retarded Hubbard U also the electron-plasmon coupling 
is stronger, as A = y^{Uoo — Uo)ujo/'2, and Uoo is an order of 
magnitude larger than Uq. In any case, if Uoo ^ Uq < ujq, the 
Lang-Firsov approach describes well the low energy proper- 
ties of the system. Here, we use the Lang-Firsov approxima- 
tion in a new and original way, as a low frequency correction 
to the dynamic factor of our Green's function ansatz. 

The factorization in the r space implies a convolution in the 
Wn space. 



G{iuJn) = -^ X! C!o{iuJm)F{iujn - iuj,n), 



(7) 



where F{ivn) are the Matsubara components of the Bose fac- 
tor. A way to improve upon the DALA is to choose the F 
factor such that the Lang-Firsov behavior is obtained at low 
frequency. This can be done, by introducing an enhanced F 



defined as: 



Fr 



,{iVn) = 



if i^„ 







1 -5exp(-^))FDALA(«i^n) elsewhere ' 

(8) 
where a, b, and c are parameters determined by the following 
conditions: F{Q) = F{/3) = 1, which keeps the correct num- 
ber of particles provided by Gstatic; G{iuJi) = GlfI^i^/) for 
I = (the first Matsubara frequency), which gives the correct 
Friedel sum rule fulfilled by the Lang-Firsov approximation 
(and broken by the DALA); c is the crossover frequency be- 
tween the LP behavior at low energy and the atomic DALA 
behavior at high frequency, and its optimal value is ?» ljq/IO, 
with cjQ the lowest Holstein frequency. 

The LF Green's function is given by the usual expression: 



GLF(iw) = 



exp 



j+Z^— A2/cjo-cxpf -^ ) A{iLo)-Yl[Uo]{iuj) 



, (9) 



which clearly implies that an Anderson model with static 
Hubbard Uo must be solved, hybridized via a renormalized 
bath exp(-A^/a;Q)A(iw). Therefore, to get the DALAh-LF 
Green's function one has to solve two static models (with reg- 
ular and renormalized bath) and mix them together by using 
our definition in Eq. |8] 

One should note that the Lang-Firsov Green's function in 
Eq. |9] has been written for a single plasmon (or phonon) 
Anderson-Holstein model, which gives a rough representation 
of the dynamically screened U{iv) present in ab-initio mod- 
els. For instance the c-RPA approximation for U{iv) usually 
leads to a quite broad spectrum of screening plasmons, which 
is difficult or impossible to fit accurately by a single frequency 
model. However, a generalization of the Eq. |9]can be easily 
done by following the same lines as in the Ref. 2lL If the LF 
renormaUzation factor exp(— A^/ojq) is replaced by 



exp 



duj Im[/(a)) 



(10) 



where U{uj) = U{u}) ~ Uq, and a) are real frequencies, all the 
screening plasmons are treated on the same footing. Although 
the LF approximation is accurate only for large uj, the Eq.fTOlis 
a good approximation if ImU (uj) / ui'^ goes rapidly to zero for 
small UJ. For the Holstein single-mode interaction lmU{u}) = 
— 7rA^((S(a) — Wo) — S{u} + ojq)), and one recovers the standard 
LF renormalization factor 

The improvement of the DALAh-LF correction with respect 
to the simple DALA factor is apparent in Fig. 13 a), to be com- 
pared with Fig. EJa). In particular, at low frequency the LF 
is capable to recover the pinning condition (G (iO+) = — 2 
for the semicircular density of states at half filling), violated 
by the DALA at low temperatures. This is reflected in the be- 
havior of the Bose factor Fdala+lf drawn in Fig.|5jb), whose 
value at /?/2 goes correctly to 1 for /? — > oo, in contrast to the 
-Fdala which does not change with temperature. Indeed, the 
exact F factor shows a quite strong temperature dependence, 
as reported in Fig.|4fd). In Fig.|3c), we plot the spectral rep- 
resentation of the DALAh-LF Green's function at /? — 160 
obtained by the DALA factor assisted analytic continuation 
described in Sec.HIl 



B. Diagrammatic first-order expansion: Gaussian cumulants 

(GO 

Another way to improve upon the DALA is to rely on the 
diagrammatic expansion of the interaction. This can be done 
in various ways. Here we took two routes: making a cumu- 
lant Green's function expansion in the full U{iiy) and the re- 
tarded part U{iv) = U{iv) — Uq only. In the former approach, 
dubbed "Gaussian cumulants" (GC) and described in this Sub- 
section, the perturbation theory is built on the Gaussian action 
containing the hybridized Qq as the bare Green's function. 
In the latter method, named "instantaneous bold cumulants" 
(IBC) and introduced in Subsection llV CI the perturbation ex- 
pansion requires the calculation of density-density correlators 
within the static Uq model by means of the CTQMC algorithm 
(or other algorithms suitable for static interacting models).' 

The first-order expansion in the full U{iv) leads to the fol- 
lowing expression for the Green's function: 



GuiT) = 00 (r) exp 



/3 



dtdt'x{t,t',T)U{t^t') 



fi — A{iuj), and xl^i^',''') 



(11) 



where Qq ^{iuj) — u 
{TcrclN{t)N(t'))Q/gQ{T) is the connected part of the 
density-density correlator easily computed for the Gaussian 
Qo propagator N{t) ~ J^a ("-ctI^) ^ 1/2) is the spin (and or- 
bital) integrated density. Note the exponential form in Eq. [TT] 
which comes from the first order cumulant expansion in U. 
According to the definition in Eq. [3j the factor F obtained in 
this way is given by the ratio Gu {t)/Guo (''')' which reads 



F,c(t) = exp -- / / dtdt'xit, t', T)U{t - t') 



(12) 



with U{t) — U{t) — UoS{t). The above expression for the 
dynamic factor has the advantage to be very accurate at both 
small U, since the interaction is treated perturb atively, and 
strong coupling, as the density-density correlator will factor 
up by giving the exact dynamic atomic limit (see the Ap- 
pendix). Incidentally, this is the reason why the cumulant ex- 
pansion is more effective than the regular perturbation theory 
for this case. 



C. Diagrammatic first-order expansion: instantaneous bold 
cumulants (IBC) 

One can go beyond the approximation for F in Eq. [12] 
and compute the density-density correlator x by using the 
interacting static Uq Green's function Go(t) as propagator 
instead of the hybridized Green's function Qo{t). This re- 
quires a correlated method, as the corresponding action is 
no longer Gaussian and the Wick theorem cannot be applied. 
The QMC methods can compute x directly in the interacting 
systems.' This development represents a consistent diagram- 
matic first-order expansion in U{iv) where the reference the- 
ory is the static C/q model. The main difference with respect 



to the previous method reported in Subsec. IIVBI stems from 
the fact that not only hybridization effects but also the im- 
pact of the instantaneous Uo6{t — t') interaction are included 
in X- The cumulant expansion of the full Green's function 
in the retarded part lJ{iv) with Go(t) taken as the instanta- 
neous bold propagator gives directly a factorized form of the 
type in Eq. [3] At the first order, it takes the expression in 
Eq. [TT]with Qq{t) replaced by Go(t) and the thermal aver- 
ages (• • • )o computed for the static interacting model, such 
that Xgo(*'^''^) = {TcrclN{t)N{t'))GjGQ{T). It turns 
out that in this case the F factor reads; 



-FibcIt) = CXp 



/3 



dMi'xG„(t,i',r)t7(i-i') 



(13) 
In the following we report a comparison of the different 
approaches proposed in this work. As common benchmark, 
we chose to compute the Green's function of the Anderson- 
Holstein impurity problem both in the symmetric and asym- 
metric cases, and check the BFA solutions against the numer- 
ically exact one provided by Rubtsov's CTQMC algorithm. 



V. OVERVIEW ON THE PERFORMANCE OF THE 
PROPOSED FACTORIZATION APPROXIMATIONS 

The factorization introduced in Eq. |3]is an extremely use- 
ful Green's function ansatz to compute thermal and spectral 
properties coming from a generic retarded interaction in the 
multiband Anderson model, and also the Hubbard model at 
the DMFT level. The product in the time domain between a 
Bose factor F embedding the dynamic properties of the in- 
teraction and an auxiliary Green's function Go including in- 
stantaneous interactions and low energy features, leads to a 
deconvolution of the spectrum into low and high frequency 
contributions. Its low frequency part, depending on Go, can 
be easily obtained by available and well developed ME meth- 
ods, while the high energy features, as the plasmon satellites, 
difficult or impossible to obtain by standard analytical con- 
tinuation, are directly given by the analytically known Bose 
factor F, which can be inverted in an accurate way by the 
Pade approximants. 

The Green's function Bose factor is a very general ansatz. 
However, the factor F is not exactly known in the generic case 
and needs some approximations. The most practical and phys- 
ically insightful one is borrowed from the dynamic atomic 
limit, where the factorization is exactly given by a Go de- 
pending only on the instantaneous Uq times a factor which 
depends only on the retarded part U{iuj). The DALA pro- 
posed in Sec.|III]consists of keeping the atomic dynamic fac- 
tor F and taking the Go from the exact numerical solution of 
the instantaneous Uq model. As reported in Figs.[T]and|6jb), 
the DALA performs well in the antiadiabatic regime when 
<^o > Uoo — Uq, and in the strong coupling regime, with a 
large Uq. The quality of the low energy part of the DALA 
Green's function worsens as the interaction becomes weaker 
and luq gets smaller, while its high energy tails are well repro- 
duced even in the intermediate coupling. The major failure of 



the DALA is the breaking of the Friedel sum rule, which is 
apparent at low temperature (below f3 — 40) or away from the 
antiadiabatic regime. That is not surprising since the DALA is 
built upon the atomic limit. In general, the DALA works when 
the high energy (unscreened) part of the retarded interaction is 
well separated (in frequency) from the low energy (screened) 
part, which is the most common situation in the case of re- 
alistic Hamiltonians, with U{iuj) determined ab-initio by the 
c-RPA approach^. 
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FIG. 6: (color online) Symmetric Anderson model with retarded 
screened interaction as in Fig. [3] ujq = 5, Uoo — Uo ~ 5, P = 10, 
for /3 = 20 and two values of Uo, (a) f/o = 1; (b) Uo = 2. The 
numerically exact Rubtsov's CTQMC, DALA, DALA+LF, GC, and 
IBC Green's functions are reported. 

An improvement upon the DALA is represented by the 
Lang-Firsov approximation, which fulfills the Friedel sum 
rule and cures the low-energy low-temperature behavior The 
way to incorporate the LP into the DALA factor, described 
in Subsec. |IV A| clearly improves the low-energy features of 
the DALA Green's function, as one can see in Figs. |6] and 
H where the value of G(^/2) yielded by the DALAh-LF ap- 
proach is closer to the exact numerical result in both the half- 
filling and quarter- filling Anderson-Holstein impurity models. 
The DALAh-LF factor provides an overall better agreement 
with the numerically exact Green's function computed by the 
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FIG. 7: (color online) Asymmetric spin unpolarized Anderson model 
with retarded screened interaction as in Fig. [6] ljo ~ 5, Ux — Uo = 
5, for 13 ^ 20 and two values of Uo, (a) Uo = 1; (b) Uo = 2. 
In this case the number of particles Na was set to be around the 
quarter filling. The numerically exact Rubtsov's CTQMC, DALA, 
DALA+LF, GC, and IBC Green's functions are reported. 



density correlator {TctCqN {t)N {t'))Q. In the first one (the 
Gaussian cumulants or "GC") the thermal average is com- 
puted with a Gaussian action based on the hybridized non- 
interacting Green's function, while in the other one (the in- 
stantaneous bold cumulants or "IBC") the brakets are com- 
puted with the instantaneous Uq Green's function. The GC 
works well in the weak interacting regime (small Uq) and it 
worsens as Uq is getting larger, as shown in Fig.|6] Away from 
half-filling, the GC Green's function is reliable in the "empty" 
part of the spectrum (for r < /3/2), and where the dynamic 
part is more relevant in setting the tails of the Green's func- 
tion. On the other hand, the IBC performs well everywhere in 
all cases when ujq is large, as seen in Figs. |6] and |7] and it is 
supposed to work quite well even away from the antiadiabatic 
regime, since it retains the feedback of the instantaneous in- 
teraction on the dynamic part via the x correlator. The price 
to pay in the latter case is the statistical uncertainty of Fibc, as 
X must be evaluated by means of Monte Carlo techniques, in 
contrast to the GC approach where x is known up to machine 
precision thanks to Wick's theorem. This could lead to some 
inaccuracy in the analytical continuation of the IBC Bose fac- 
tor. 

If the goal is to compute the spectral representation of a 
Green's function with retarded interactions, we found that 
the most effective way is to use the Bose factor taken from 
the dynamic atomic limit. As already reported in Sec. HI] 
one has to define an auxiliary Green's function G^ux{t) as 
G{t)/Fuala{t), and then use Eq.|4]to compute the spectrum. 
The DALA factor is the simplest to invert, and its bosonic 
spectral representation is the "physical" density of plasmonic 
(or phononic) modes, whereas the other factors beyond the 
DALA include also some low energy contributions. There- 
fore, the DALA factor is the most recommended in the as- 
sisted analytical continuation. 



VI. APPLICATION TO THE SINGLE-BAND 
HUBBARD-HOLSTEIN MODEL 



CTQMC method. The LP correction turns out to be important 
particularly in the weak-intermediate coupling and at low tem- 
perature, while in other cases the simple DALA approxima- 
tion is already good enough. For instance, in Fig. |6]b), when 
Uq is strong and the temperature not-so-low (/3 = 10), the 
DALA Green's function coincides with the DALAh-LF one, 
and both are on the top of the exact numerical solution. A 
common limitation of the DALA and DALAh-LF approxima- 
tions is that they become inaccurate when ujq < Uoo — Uq, 
while the work very well in the antiadiabatic regime. 

A different route to determine F is built on diagrammatic 
expansion techniques. We have proven that the cumulant ex- 
pansion of the Green's function up to the first order in the in- 
teraction provides the exact connection to the dynamic atomic 
limit (see Appendix). Indeed, the F factor computed in this 
way (see Eq. [T2] | has the exact atomic behavior in the zero 
hybridization limit. We propose two flavors for the cumu- 
lant expansion. In both cases one has to evaluate the density- 



In the following, we are going to present applications where 
our BFA approach is used as solver of the Anderson impu- 
rity problem resulting from the DMFT self-consistency con- 
ditions. Therefore, in contrast to what has been shown in the 
methodological Sections, the following results are the con- 
verged solutions of the DMFT equations for the full lattice 
Hamiltonian, which incorporates the retarded f/(w) as the on- 
site interaction. 

The first application is on the single band half-filled 
Hubbard-Holstein model solved by the DMFT on the Bethe 
lattice. To study the impact of the retarded screened interac- 
tion on the Mott transition, we analyze the half-filled model at 
/3 = 40 with different Uq, Uoo, and ujq parameters. We choose 
to work always with Uoo > Uq, a "physical" condition which 
states that the unscreened U is larger than the screened one. 
The screening frequency is taken such that we are in the antia- 
diabatic regime, and we use our various Bose factor approxi- 
mations to predict the critical value of Uoo for the Mott tran- 
sition, once the other parameters are fixed, such that luq — 10 



and Uq = 2. As reported in literaturer^ii^ the metal-Mott 
insulator transition in the Hubbard-Holstein model is first or- 
der. In Tab. I] we present our lower and upper critical values 
obtained for /? = 40 by means of the BFA at various levels 
of approximation. As one can see, the actual values at criti- 



TABLE I: Upper and lower critical Uoo values obtained at /3 = 40 
by different BFA methods for the Mott transition in the single-band 
Hubbard-Holstein model on the Bethe lattice at half-filling, once 
Uo{= 2) and ijjo{= 10) are fixed. The values are expressed in the 
half-bandwidth units. By our BFA we always get a first order phase 
transition, with a hysteresis effect, in agreement with what reported 
in literature'*. 
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cality depend quite strongly on the approximation used, with 
the DALAh-LF giving the results aligned with the numerically 
exact CTQMC method by Werner and Millis^' (see Fig. |8]l. 
In order to have a close comparison to the data published in 
Ref. |2ll, we caiTied out some calculations also for /? — 100. 
The agreement between the DALAh-LF and the Monte Carlo 
predictions is quite remarkable. This highlights again the im- 
portance of the LF coiTection at low energy in order to accu- 
rately predict the physical properties. 

Fig. [8] summarizes the results for the upper critical line for 
different sets of ojq and Uq, whose values are taken not so 
far from the critical Uc2 ~ 2.6 of the static model computed 
at /3 = 4022i^. It is clear that the fully retarded model with 
Uoo > Uq is more correlated than the static one with the same 
instantaneous Uq. The Mott transition happens at values of 
screened Uq lower than the critical static 11^2 for any finite luq. 
The dependence of the critical parameters on ujq is also clear 
At fixed screened Uq, a smaller ujq corresponds to a smaller 
unscreened Uoo at which the Mott transition is reached. In- 
deed, if the frequency of the screening plasmon is closer to 
the Fermi level, it is easier for the unscreened part to induce 
the Mott transition at low energy. The same conclusions were 
reached by Werner and Millis by means of their numerically 
exact CTQMC algorithm. This shows the impact of the dy- 
namic screening features on the low energy properties of the 
model. In order to get the same effective low energy parame- 
ters, the "effective" instantaneous Uq can be up to 20 — 25% 
larger than the true screened value even for plasmon frequen- 
cies ujq{^ Uq,^ D) in the antiadiabatic regime. 

In Fig.|9]we report the spectral representation of a half-filled 
Holstein-Hubbard model with ijJq — \Q,Uq — 2, Uoo = 6.5, 
and temperature /3 = 100, quite close to the Mott transition. 
The spectral function has been obtained with the help of the 
DALA factor by the method in Sec. [Ill It describes with high 
accuracy not only the low-energy features but also the plas- 
mon satellites in the anti-adiabatic regime, typical of realistic 
materials. 
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FIG. 8: (color online) Critical Ud line of the Mott transition calcu- 
lated at /3 = 40 and half-filling on a Bethe lattice for the screened in- 
stantaneous interaction f/o as a function of Uoo/Uo and ljo ~ 10, 20 
(green diamonds and red circles respectively). The critical values 
are obtained from the Green's function computed by the DALA+LF 
ansatz. The behavior of the critical screened U is in agreement with 
the one computed by Werner and Millis with their CTQMC algo- 
rithm at LJo ~ 10 and /3 — 100 (black upper triangles). Indeed, the 
critical line calculated by means of our DALA+LF approach (blue 
lower triangles) is on top of their CTQMC points"'. This highlights 
the accuracy of our approximated DALA+LF method if compared to 
the exact numerical result. Note the quite strong temperature depen- 
dence of the critical U^i in going from P — 40 to j3 — 100. This is a 
quite interesting effect which deserves further analysis. 
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FIG. 9: Spectral representation of the Green's function for a 
Holstein-Hubbard model with loq — 10, Uq = 2, Uoo ~ 6.5 and 
/? — 100, close to the blue Mott transition line of Fig. [8] The 
analytic continuation has been assisted by the DALA Bose factor, 
which allows one to describe accuratley the plasmon satellites cen- 
tered around frequences multiple of luq. 
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VII. REALISTIC APPLICATION TO SrVOa 

Our second application is on a realistic Hamiltonian for 
SrVOs, a very well studied material which represents a bench- 
mark for theories describing strongly correlated compounds. 
SrVOs is the prototype of a correlated metal where the many- 
body treatment of correlation in the d manyfold is important 
to explain the spectral properties^. Therefore, it has been the 
subject of intensive studiesi^lri^^ applying DMFT in the con- 
text of realistic strongly correlated Hamiltonians. Indeed, its 
band structure is relatively simple due to its undistorted per- 
ovskite structure, resulting in the occupation of one electron 
in three degenerate t2g bands crossing the Fermi level. The 
p Oxygen ligands are quite well separated from the d levels, 
such that the definition of a low-energy 2g-Hamiltonian is un- 
ambiguous. Thus, SrVOa has been the testing case for many 
new DFT-DMFT implementations. ^'^^^^ On the other hand, 
SrVOa has been the subject of intensive experimental activity, 
with magnetic, electrical, optical measurements^^^, and by 
means of photoemission spectroscopy (PES)^^ and angled 
resolved PES (ARPES)'*"^ 

Here, we consider a model where only the t2g electron is 
retained, and all the others contribute to screen the local lattice 
interaction. The DFT band structure has been calculated with 
the linear muffin-tin orbital (LMTO) framework in the atomic 
sphere approximation (ASA), which allows one to work with 
a native d-projected and localized orbital representation. The 
realistic retarded U for this compound has been computed in 
Ref.|5| based on the c-RPA construction. 

The low-energy Hamiltonian we are going to work with 
consists of the LDA 2g-Hamiltonian, the plasmonic part giv- 
ing rise to the dynamical screening and the following static 
interaction Hamiltonian: 



U' 
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where nm(y=d\y^„d^n„ is the usual density operator, with m,(j 
denoting orbital and spin indexes, and d\^^, d^^ represent- 
ing the i2g localized orbitals. For the 2g-orbitals of SrVOa, 
U' = U -2J a.ndU" = f/-3J, with the screened value of the 
interaction Uo = 3.6 eV, and the Hund's coupling J = 0.68 
eV. The additional retarded interaction which is included in 
our model, couples to the total charge of the system, as de- 
scribed by Eq. |2] There is no need for an explicit double- 
counting term, since such a correction is absorbed into the 
effective chemical potential fixing the particle number to one. 
To study the impact of C/retarded on the low energy prop- 
erties of the model, we took into account different Holstein 
single-plasmon C/(ii^)'s, and compared them with the corre- 
sponding static model. We used our BFA approach in its 
simplest DALA formulation, and computed the spectrum at 
/3 = 10 and /3 = 20 eV"^ Indeed, it turns out that at 
those temperatures the LF correction is irrelevant, and so the 
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FIG. 10: (color online) Spectral function of one electron in 3 i2g 
bands for the SrVOa obtained by the DFT in the local density ap- 
proximation. The DMFT calculations are done a.i j3 — 10 eV~^ 
with the DALA Bose factor solver for the retarded U and the Hirsch- 
Fye QMC solver for the instantaneous interaction. Three types of 
interactions are reported: the purely static one, with Uo = 3.6 eV, 
a dynamic Holstein one with the same instantaneous interaction but 
with the unscreened Uao = 7 eV and a plasmon frequency a;o = 5 
eV, another Holstein interaction with the same Uo and Uoo but dif- 
ferent u}o{= 15 eV). Note the displacement of the upper Hubbard 
shoulder shifted at lower energy by the presence of the plasmon at 
cjo = 5 eV. 



DALA performs here at best even in terms of efficiency. In 
Fig. [To] we report the spectra computed with the instanta- 
neous Uq{— 3.6 eV) interaction, and two retarded interac- 
tions with the same unscreened Uoo{— 7eV), but different 
screening wo's (one at 5 eV, the other at 15 eV). The static 
model has a quasiparticle peak at the Fermi level typical of 
strongly correlated compounds, with a lower and upper shoul- 
ders reminiscent of the lower and upper Hubbard bands. In the 
dynamic model, the effect of the 15 eV plasmon is to renor- 
malize the quasiparticle spectral weight, transferred at higher 
energies, while the shape of the low energy spectrum is almost 
unchanged. The effect of the 5 eV plasmon, closer to the low 
energy sector, is more remarkable. Beside a stronger spectral 
weight reduction, there is a shift of the upper Hubbard shoul- 
der to lower energies than in the corresponding static model. 
This effect is certainly due to the interplay between the plas- 
mon satellite at 5 eV, visible in Fig. [TOl and the low energy 
features of the spectrum. 

The realistic retarded U{iv) is characterized by two main 
screening frequencies at 5 and 15 eV^'', while the unscreened 
value of C/ is 16 eV, much larger than the one in the mod- 
els considered so far. By applying our Bose factor ansatz 
to this problem, we found the spectral function in Fig. [TTl 
reported together with the instantaneous interaction with the 
same screened C/o(= 3.6 eV). Also here we note the quasi- 
particle weight reduction, and the shift of the upper Hubbard 
shoulder to lower frequencies, as in the model U analyzed 
before, while the lower Hubbard band is almost unchanged 
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FIG. 1 1 : (color online) Spectral function at /? = 10 eV~ of one 
electron in 3 t2g bands for the SrVOs obtained by the DFT in the 
local density approximation. The DMFT calculations are carried out 
with the DALA Bose factor solver for the retarded U and the Hirsch- 
Fye QMC solver for the instantaneous interaction. The result of the 
realistic Green's function is reported, showing a more correlated be- 
havior than the corresponding static model with the same instanta- 
neous Uo{= 3.6 eV). The quasiparticle peak is smaller than the 
corresponding static model, with a spectral weight transfer at higher 
energies. 



with respect of the static model. Its maximum corresponds 
to the position found in PES measurements, and well doc- 
umented in previous DMFT studies on SrVOs. The weight 
reduction of the quasiparticle peak coincides with a smaller 

value of Z — 1/ (l — dT,' {iuj)/duj ), namely a larger 

value of the effective mass m* . In particular, for the realistic 
dynamic model we found a. Z k, 0.5, which gives an effective 
mass renormaUzed by 2 with respect to the DFT band struc- 
ture, while for the corresponding static model we obtained a 
value of Z « 0.7. Recent ARPES data yielded an effective 
mass m* ss 2to[ |"^^'^^'^' , which is in a good agreement with 
our findings for the realistic retarded interaction. On the other 
hand, the static model with the same screened C/q underes- 
timates the correlation by a factor of 1.4, and so the value 
of m*. In the literature, a static model with a larger instan- 
taneous C/o(« 5 eV) has been used to find the experimental 
mass renormalization^^. Such a larger value of U could be 
justified by the constrained LDA method''^-^'* used to deter- 
mine a priori the on-site interaction, but known to overesti- 
mate its strength. The difficulty here is to reproduce by the 
same model both the effective mass and the position of the 
lower Hubbard band, which turns out to be shifted at lower 
energies (w — 2 eV) by a stronger C/q- This is the reason why 
in Ref. 37 the authors made the choice to work with C/q = 4 
eV, a slightly weaker effective static interaction which gives 
the lower Hubbard band correctly peaked at —1.5 eV. Some 
cluster calculations give the same position for the lower Hub- 
bard band, with the interpretation that its peak depends on 
the hybridization and screening properties provided by the 
ligand s^^-^^ . With our dynamically screened model we de- 



scribe correctly both the mass renormalization (2mo) and the 
lower Hubbard band position peaked at —1.5 eV, as apparent 
in Fig. [TT| This analysis highlights the importance of includ- 
ing the proper retarded interaction, to have a reliable and fully 
ab-initio description of the correlation in these materials. 

We note however, that the present description includes only 
the 2g orbitals, so that the above conclusions are valid in an 
energy range where no other orbital contributions are present. 



VIII. CONCLUSIONS 

We introduced a factorized form G — FGq for the Green's 
function of the Anderson model with generic retarded (dy- 
namic) interaction, dubbed "Bose factor ansatz". We pro- 
posed various approximations for the Green's function Bose 
factor F, the most practical and effective one borrowed di- 
rectly from the dynamic atomic limit (DALA), whose form is 
analytically known. The DALA provides an accurate way to 
compute the Green's function in the antiadiabatic limit and 
evaluate the spectral properties in the full frequency range 
by means of an improved analytical continuation method. In 
practice, the inversion from the imaginary to the real fre- 
quency domain is assisted by the factor F which retains the 
main information on the position and strength of the plasmon 
satellites, and enters in the analytical continuation as a con- 
volution with the instantaneous part Gq. We carefully an- 
alyzed the pros and cons of the DALA, and found various 
ways to improve F, by either using the Lang-Firsov approx- 
imation at low-frequency, or resorting to diagrammatic tech- 
niques. Finally, we applied our approach to lattice problems 
in the context of the DMFT formalism. We took into account 
the Hubbard-Holstein Hamiltonian at half-filling and on the 
Bethe lattice, and we studied the Mott transition in the spirit 
of looking at the retarded interaction as resulting from "real- 
istic screening" of the bare [/ by a single plasmon. The sec- 
ond application has been for the realistic SrVOs Hamiltonian, 
where the t2g electrons interact via a retarded on-site U previ- 
ously determined ab-initio at the c-RPA level'*. In both cases, 
it turns out that it is important to retain the retarded features 
of the local interaction resulting from the dynamic screen- 
ing in order to have a reliable ab-intio description of corre- 
lated materials. Also, our approach could be useful to deter- 
mine whether some spectral signatures at intermediate energy 
(« 10 eV) seen in a broad class of correlated materials^^ come 
from a bulk dynamic screening. In perspective, more work has 
to be done theoretically to rationalize the effects of the screen- 
ing on both the ground and excited states properties of corre- 
lated compounds. Morever, by means of the same formalism, 
one can study Jan-Teller models to describe the impact of the 
electron-phonon coupling to the spectral function of the dis- 
torted compound. Last but not least, dealing with a frequency 
dependent U is an essential step toward the implementation 
of the GWh-DMFT framework^^, where the screening result- 
ing from the GW polarization has to be included consistently 
in the low-energy correlated model solved at the DMFT level. 
Therefore, behind this work there are important experimental, 
theoretical and methodological implications that one can now 
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start taking into account. 



Appendix: From the first-order cumulant expansion to the 
dynamic atomic hmit 

In this appendix we prove that the Bose factor F reported 
in Eq. [12] fulfills the exact atomic limit, i.e. it equals the one 
in Eq. |6]for A = 0. This is a non trivial property, which guar- 
antees the GC approximation to provide a reasonable descrip- 
tion of the insulating phase (and the metal-to-insulator transi- 
tion) in the DMFT framework, where the hybridization func- 
tion of the Anderson impurity assumes a crucial frequency 
dependence in order to represent the coupling with the self- 
consistent bath. In the strong coupling case, the hybridization 
A{iu!n) goes to zero for small a;„. Thus, in the DMFT lan- 
guage, this phase is mapped into an Anderson impurity prob- 
lem close to the atomic limit, and its accurate solution is re- 
quired around A = 0. 

We start by noting that the density-density correlator x fac- 
torizes in the atomic limit, as 



xit,t',T)^^x\t,r)xHt', 



(A.l) 



where x^(i,T) — (7'cT-cjA^(i))o/^o(''')- In the most general 
case, X can be resolved into its spin and orbital components, 
by defining x<r<.'(i,i',r) = {rcrclN,{t)N,,{t'))o/go{T), 
with Na-{t) = n„{t) — 1/2. The factorization in Eq. lA.ll holds 
also for the spin resolved quantities: Xaa' = XaXa' ■ Once the 
Wick theorem is applied and the former correlators are writ- 
ten in terms of the Qo{Tys, it is straightforward to prove the 
spin resolved identity and consequently the integrated one of 
Eq. lA.ll by using the atomic limit expression for Go{t). For 
instance, in the atomic limit and for a generic fi, the connected 
part of x^ reads: 



xHt, 



-e-i'P/ (1 + e 
l/(l + e-^^) 



-M/3^ 



for T < t < ^ 
for < t < T 



(A.2) 



Now we use a Hubbard-Stratonovich (HS) transformation 
to rewrite the Fee factor in a form which includes x^ in linear 
terms only. We obtain the following identity: 



where J Vcj) is the functional integral over the complex HS 
field (j>{t). By exploiting the atomic limit expression of x^ in 
Eq. IA.21 one can compute the integral involving the product 
of ip and x^, which gives the result: 



dtx\t,T)(b{t) 



'fJ.P 



1 + e-t^f^ 






1) , (A.4) 



with (j)n — /q dt(j){t)e^'^"^ the Fourier components of the HS 
field. By integrating the functional integral in Eq. IA.3l in the 
Fourier space, one gets the final expression for the exponent 
of the GC factor: 



/3 



dtdt'x{t,t',T)U{t~t') = 



-l^PtJiiv^) 



-t^l3 



-y 



1 + e-t'P 



P 



A=0 
2 



(A.5) 



We note here that U vanishes in the static limit, as U{ivn) — 
U{ivn) — Uq. Therefore, the first term in the right-hand side 
of the above equation vanishes as well. We are left with the 
last term in Eq. IA.5l which is exactly equal to the exponent of 
the DALA factor. Thus, we have proven that 



i^Gc(r)^=^FDALA(T), 



(A.6) 



for an arbitrary /i and U. 

To conclude, the fact that the first-order cumulant expan- 
sion fulfills the exact dynamic atomic limit justifies the use of 
cumulants instead of the standard first-order developements, 
and validates also the cumulant of the instantaneous bold fac- 
tor in Eq. [T3] 



1 



/3 



exp|--yy dtdt'xHt,r)U{t-t')x'{t\T)\ = 
P0exp [ -i / dtx^(i,T)0(t) 

/3 \ 



dtdt'(t){t)U-\t - t')(j){t') , (A.3) 
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